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(57) Abstract 

A method and system of detecting mass to charge ratio of ions. The method includes producing charged ions in a vacuum, accelerating 
the charged ions in an electric field into a free flight tube and detecting the charged ions at a detector associated with the free, flight tube. 
A control system selects a bandwidth for filtering a signal produced by the detector and the signal produced by the detector is then filtered 
with a variable width digital filter based upon the selected bandwidth. The bandwidth for filtering the signal may be selected from a look-up 
table within the control system based upon the mass to charge ratio of an ion of interest. Alternatively, a peak bandwidth within the signal 
produced by the detector may be determined and the signal produced by the detector may then be filtered with the variable width digital 
filter based upon the determined peak bandwidth. 
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VARIABLE WIDTH DIGITAL FILTER FOR TIME-OF-FLIGHT 

MASS SPECTROMETRY 

5 This application claims priority from U.S. Provisional Patent Application 

Serial No. 60/134,072 (Atty. Docket No. 01 6866-003400), filed May 13, 1999, the 
disclosure of which is incorporated herein by reference in its entirety. 

BACKGROUND OF THE INVENTION 
10 1 . Field Of The Invention 

The present invention relates to a software filter for use with time-of-flight 
mass spectrometry, and more particularly, to a variable width digital filter for use with 
time-of-flight mass spectrometry. 
2. Description Of The Prior Art 
15 With reference to Figure 1, time-of-flight mass analyzers or spectrometers 

consist of a source/extraction region 10, a drift region 1 1 and a detector 12. In the source 
region, an electrical field (E=V/s) accelerates the ions to a constant energy. The drift 
region is field free and ions cross the drift region with velocities that are inversely 
proportional to the square root of their respective masses. Thus, lighter ions have higher 
20 velocities and arrive at the detector sooner than heavier ions. 

In the ideal situation where ions are formed at a single point in the source 
region, ions are accelerated to the same final kinetic energy: 

1 ^ v 
— mv = eV 

2 



25 



and cross the drift regions with velocities: 

lev 

v - 



m 

and flight times: 

-|l/2 

m 

t = 

2eV 



These relationships depend upon the square root of the ions' masses. 
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In a mass spectrometer, the mass resolution is defined as m/Am. In a time- 

of-flight mass spectrometer in which ions are accelerated to constant energy: 

m _ t 
Am 2&t 

In time-of-flight mass spectrometers, it is not unusual to see a wide mass 
5 range being scanned at any given time. Ions with molecular weights between 100 and 
several thousand Da, ions ranging from 3,000 to about 20,000 Da, as well as all ions 
greater than 20,000 Da, are typically simultaneously studied in such techniques as Surface 
Enhanced Laser Desorption Ionization SELDI and Matrix Assisted Laser Desorption 
Ionization (MALDI). 

10 The fundamental physical processes involved in the previously mentioned 

processes are such that signals created by heavy ion populations are generally composed 
of lower frequency components than their light ion counter parts. For signals created by 
light ion populations, broad detection bandwidths are required to accurately sample these 
fast transients allowing for enhanced resolution mass measurement. . Signals from heavier 

15 ion populations typically do not possess significant high frequency components and thus 
may be sampled at significantly lower bandwidth frequencies. Table 1 lists theoretical 
major frequency components and estimates peak widths and mass resolutions of various 

ion signal populations along with their estimated times of flight and molecular weights as 

generated by a SELDI or MALDI time-of-flight mass spectrometer with one-meter drift 

20 region and 25 keV total energy. 



TABLE 1 



ion Molecular 
Weight (m/z) 


Ion Flight Time 
(uSec) 


Major Component 
Frequency (MHz) 


Peak Width At Half 
Height (uSec) 


Mass 
Resolution 


500 


10.2 


740 


0.0010 


5000 


1,000 


14.4 


500 


0.0016 


4500 . 


2,000 


20.4 


250 


0.0034 


3000 


5,000 


32.2 


70 


0.0134 


1200 


15,000 


55.8 


19 


0.0254 


1100 


40,000 


91.1 


2 


0.3037 


150 


150,000 


176.3 


.290 


1.7600 


50 


250,000 


227.6 


.130 


3.8000 


30 


500,000 


321.9 


.063 


8.0500 


20 
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Thus, by reviewing Table 1, it can be seen that the peak width at half 
height and mass resolution of a given ion population can be correlated to ion flight time 
for a given ion total kinetic energy and a given free flight distance. Most time-of-flight 
mass spectrometers incorporate a fixed drift region distance. Furthermore, these devices 
5 also operate using either a fixed level or precisely selectable levels of ion acceleration, 
thus allowing qualified approximations of ion total kinetic energy. Under such 
conditions, it would be possible to predict the signal frequency requirements for a variety 
of ion populations based upon their time of detection. 

The wider a peak width is, the more ions of different "sizes" may be 
10 contained within the particular ion population that is being detected. Hence, it is 
desirable to accurately display peak widths. 

Just as in other forms of spectroscopy, time-of-flight mass spectrometry 
has several sources of signal noise. Such signal noise may increase peak widths. Typical 
noise sources such as sampling noise (alaising), Johnson noise, and flicker noise 
15 contribute to the total system noise. However, sensible engineering approaches will often 
reduce these noise sources to insignificant levels. Often, the most frequently encountered 
noise in time-of-flight mass spectrometry measurements is high frequency noise created 
by the detection apparatus. The combined u se of secondary ions/electron pero ration 
schemes with high gain electroemissive detection surfaces frequently introduce high 
frequency noise that is the direct result of spurious background gas ionization, thermal or 
low energy photon noise (dark current noise), as well as higher energy photon or other 
particle-induced noise. Thus, when considering the above factors regarding ion signal 
component frequencies and time-of-flight mass spectrometry noise characteristics, it is 
evident that a fixed width filter is not a desirable solution for addressing noise problems. 
A filter in which the bandwidth may be varied over time range of the time-of-flight 
spectrum may better optimize the tradeoffs between increasing the signal to noise ratio 
while having the least negative effect on the mass resolution. 



20 



25 



SUMMARY OF THE INVENTION 
30 A method of detecting mass to charge ratio of ions in accordance with the 

present invention includes producing charged ions in a vacuum, accelerating the charged 
ions with an electric field into a free flight tube, and detecting the charged ions at a 
detector associated with the free flight tube. With a control system, a bandwidth for 
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filtering a signal produced by the detector is selected. The signal produced by the 
detector is then filtered with a variable width digital filter based upon the selected 
bandwidth. 

In accordance with one aspect of the present invention, the bandwidth for 
5 filtering the signal is selected from a look-up table within the control system based upon 
the mass to charge ratio of an ion of interest. 

In accordance with a further aspect of the present invention, the method of 
detecting mass to charge ratio of ions further includes determining a peak bandwidth 
within the signal and filtering the signal produced by the detector with the variable width 
10 digital filter based upon the determined peak bandwidth. 

Accordingly, the present invention provides a system and method, 
especially well suited for time-of-flight mass spectrometry wherein the width of a digital 
filter of varied over the mass spectrum to optimize the signal to noise improvement 
throughout the mass range. This is done without significantly compromising the mass 
15 resolution. 

Other features and advantages of the present invention will be understood 
upon reading and understanding the detailed description of the preferred exemplary 
embodiments found hereinbelow, in conjunction with reference to the drawings, in which 



20 



like numerals represent like elements. 



BRIEF DESCRIPTION OF THE DRAWINGS 
Figure 1 is a schematic illustration of a time-of-flight mass spectrometer; 
Figure 2 is a schematic diagram of one possible embodiment of a mass 
spectrometer system in accordance with the present invention; and 
25 Figures 3-6 are graphs illustrating the effect of a variable width digital 

filter in accordance with the present invention on signals from a mass spectrometer. 



DETAILED DESCRIPTION OF THE PREFERRED EXEMPLARY EMBODIMENTS 
Generally, a mass spectrometer 20 charges or ionizes molecules of a 
30 sample S into ions P in a vacuum 21. These ions are accelerated by an electric field 

produced by an ion-optic assembly 22 into a free flight tube 23. The velocity at which the 
ions may be accelerated is proportional to the square root of the accelerating potential, the 
square root of the charge of the molecule and inversely proportional to the square root of 
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the mass of the molecule. The charged ions travel, i.e., "drift" down the time-of- flight 
tube to a detector 24. 

Detector 24 generates a signal, generally an electronic signal, that is 
generally received, and preferably stored, in a control system 25, such as for example, a 
5 computer or the like. The signal is then displayed on some type of a display screen 26, 
such as a computer monitor, an oscilloscope, etc. Such viewing may be done either in 
real time, i.e., as the signal is received from the detector, or from the stored signal. 

As previously discussed, time-of-flight mass spectrometry has several 
sources of signal noise that include, for example, sampling noise (alaising), Johnson 
10 noise, and flicker noise. Some of the most frequently encountered noise in time-of-flight 
spectrometry measurements is high frequency noise created by the detector. Often the 
detector includes secondary ion/electron generation schemes with high gain 
electroemissive detection surfaces that frequently introduce high frequency noise that is 
the direct result of spurious background gas ionization, thermal or low energy photon 
15 noise (dark current noise), as well as higher energy photon or other particle-induced 
noise. 

A digital filter is a linear/shift-variance system for computing a discreet 
output seq uence form a discreet input sequence. The digital filter is applied to a series of 
equally spaced data values f k s / (/,), where t x = t 0 + A for some constant sample spacing 
20 A and i =... -2,-1,0,1,2... The simplest type of digital filter, commonly referred to as a 
fixed width moving average filter, replaces each value /,• by a linear combination gj of 
itself and some number of nearby neighbors, 



Here ni and rtn are the numbers of data points used. to the left and to the 
25 right of data point / respectively. n L and n R are both constants, thus, the filter has a 
constant width of n L + 1 . 

Replacing constants and hr with nu = n L and nm = (r,-) creates a 
variable width digital filter. 



WO 00/70648 PCT/USOO/13153 

6 

Such a variable width digital filter 27 is included with control system 25 
and is applied to signal data received from the detector to increase the signal to noise ratio 
of the spectrum. The variable width digital filter utilizes the fact that data is over- 
sampled in the time domain. Increasing the filter width decreases the signal bandwidth, 
5 and can improve the signal to noise ratio if the signal of interest is a far lower frequency 
than the noise. Thus, the variable width digital filter offers the ability to filter the data 
after the data has been recorded without permanently changing the raw spectrum data. 
Therefore, the data may be examined and the filter adjusted to optimize the trade-offs 
between signal to noise enhancement and resolution that filtering the data imposes. 

10 In a preferred embodiment, a look-up table 28 is provided within the 

control system for selecting the bandwidth of the variable width digital filter. The 
bandwidth is selected based upon observed or theoretical data related to varying mass to 
charge ratios of ions. Alternatively, the control system may be programmed such that 
upon receiving signals from detector 24 or upon analyzing saved or recorded signals 

15 already received from the detector, the widths of the various peaks within the signal are 
determined in the bandwidth for the variable width digital filter as used therein. 

Figures 3-6 illustrate the improvement in signal to noise and resolution 

that a variable width filter in accordance with the present invention provides. The same 

spectrum is compared unfiltered, with a 51 point fixed moving average filter, and with a 

20 variable width moving average filter. The data was acquired at 250 MHz, making the 51 
point moving average filter .024 microseconds wide. The variable width filter varies its 
width in points by interpellating a table of Mz to expected peak widths. Table 2 
illustrates an example of width values that were used to calculate the variable filter widths 
for the figures. In the figures, the X-axis represents the mass/charge ratio (M/Z) while the 

25 Y-axis represents arbitrary ion intensity units. 



TABLE 2 
Variable Width Filter Table 



Position (Daltons) 


Width (Daltons) 


0 


10 


33000 


650 


147000 


5500 
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Figure 3 illustrates the entire spectra of example data. Figures 4-6 
illustrate the detailed view of peak 1, which occurred at 6,634 daltons, a detailed view of 
peak 2, which occurred at 18,123 daltons, and a detailed view of the third peak, which 
occurred at 70,567 daltons, respectively. Tables 3 and 4 illustrate the effect of filtering on 
M/Z resolution and the effect of filtering on the signal to noise ratio. 



TABLE 3 
Effect Of Filtering On M/Z Resolution 



M/Z Resolution 


Peak # 


M/Z 


No Filter 


51 Point Fixed 


Variable Width 


1 


6634 


147 


71 


135 


2 


18123 


82 


42 


63 


3 


70578 


40 


33 


42 

C3 



TABLE 4 

Effect Of Filtering On Signal To Noise Ratio 



Signal To Noise 


Peak# 


M/Z 


No Filter 


51 Point Fixed 


Variable Width 


1 


6634 


10 


157 




2 


18123 


21 


616 


424 


3 


70578 


.79 


24 


106 



Thus, it can be seen that with a fixed width filter, the optimized filter value 
for mass resolution and signal to noise enhancement can only occur at a relatively small 
portion of the spectrum. With a variable width filter in accordance with the present 
invention, the entire spectrum may be optimized. Accordingly, it is easier to isolate peaks 
and therefore isolate ions as opposed to groups of ions. 

Appendix A contains source code that provides an example of a variable 
width digital filter for time-of-flight mass spectrometry in accordance with the present 
invention written in C++. 

Thus, the present invention provides a digital filter for time-of-flight mass 
spectrometry that varies the width of the filter over the mass spectrum to optimize the 
signal to noise improvement throughout the mass range without significantly 
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compromising the mass resolution. This is accomplished by predicting the required filter 
width at a given time in the spectrum The predicted widths may be generated from 
theoretical or observed spectra. 

Although the invention has been described with reference to specific 
exemplary embodiments, it will be appreciated that it is intended to cover all 
modifications and equivalents within the scope of the appended claims. 
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APPENDIX A 



Example Source Code for a Variable WJdth Filter in C++: 

void CCFi IcerVarWidth : :FilterDaca (const CArray *doubl e , doubl e > trawData, 

C Array <double . doubl e> & f i 1 ceredDaca . 
CRangt;<long> iusabl eAange ) const 



.1 



long n * rawDa ca . GecSi ze t ) ; 
f il ceredDaca . SecSiie (n) ; 

if in «- 0) 
( 

usableRange , Setup < 0. 0) ; 
I'eturn ; 

) 

ASSERT (m^spectrum t» MULL) ; 

// Prepare an array' of filter widths. 
CArx ay <doubl e . doubl e> idenc ; 
CArray<doubl e , doubl e > width; 
CArray <1 ong , long* halfWidth; 
idenc . SetSi ze in) ; 

for (long i = 0,-i<n;+ + i) » 
idenc I i ] * i ; 

ro_varWidth->EvaluateOnCrid (m_spectrum. CCSpecFun : : IUOP . CCSpecl"un : : DUDPWid . ident. width) 

hal f Width. SetS ize (n) ; 
for (i * 0; i < n; ++iV 

halfWidthti) = (long ) (width I i J > / 2; 

// Do Che filtering. 

long haveFrom ■ 0, haveTo * -li 

double sum - 0; 

long usableFrom = 0, usableTo = n - 1; 

for (i - Q; i < n; * + i) 

-( ~ 

long hw = hal f Width ( i ) ; 
if (hw < 0) 

hw = 0; 
long from = i - hw; 
long co = i ♦ hw; 

i f ( from < 0) 
( 

from * 0; 

if (usableFrom < i t 1) 

usableFrom « i * 1 ; 

) 

it (co >= n) 
( 

co » n - 1 ; 

if (usableTo > i - 1) 

usableTo * i - 1; 

) 



if I from < haveFrom) 
( 

for (long k = from; k < haveFrom; * + k) 

sum rawData(k| ; 
haveFrom * from; 

) 

else if (from > haveFrom) 
( 

for (long k * haveFrom; k < from; + + k] 

sum - = rawData(k); 
haveFrom = from; 

) 

if (co > haveTo) 
( 

for (long k * co; k > haveTo; --k) 
» n 
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sum rawData(k) ; 
haveTo » CO; 

) 

el se i f ( co < haveTo) 

for (long k - haveTb; k > to; --k) 

sum -= rawDaca (k) ; 
haveTo = to; 

) 

f il teredData I i ) 3 sum / (to - from + 1); 

} 

// Check for the 'nothing usable' case, 
if (usableFrom >~ usableTo) 

t 

inc k ~ n/2 - K 

usableRange . low = usableRange .high * k; 
double v = f ilceredData (k) ; 
for (i = 0; i < usableFrom; 

t ilteredDataU ] - v; 

) 

else 

usableRange . low « usableFrom; 
usableRange . high = usableTo + 1; 



// Fill the tails with constant. £i 1 teredData ( 0 ... usableRange . low- 1 ) = 

£n " re ^ aCa,OS ^ 1 ircer^:la"s a M e Ran g e.h ig h...N-l) . f U t eredDa t a,usaMeRan g e. h i g h-l I 

for (i - 0; i .< usableFrom; + + i) 

f ilceredData(i ) * fi 1 teredData ( usableFrom) ; 
for (i • usableTo +1; i < n; ++i) - 

filteredDatat i ) = fi 1 teredData I usabl eTo } ; 
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WHAT IS CLAIMED IS: 

1 LA method of detecting mass to charge ratio of ions, the method 

2 comprising: 

3 producing charged ions in a vacuum; 

4 accelerating the charged ions with an electric field into a free-flight tube; 

5 detecting the charged ions at a detector associated with the free-flight tube; 

6 selecting, with a control system, a bandwidth for filtering a signal 

7 produced by the detector; and 

8 filtering the signal produced by the detector with a variable width digital 

9 filter based upon the selected bandwidth. 

1 2. The method of claim . 1 wherein the bandwidth for filtering the 

2 signal is selected from a look-up table within the control system based upon the mass to 

3 charge ratio of an ion of interest. 

1 3. The method of claim 1 further comprising: 

2 determining a peak bandwidth within the signal produced by the detector; 

3 and ' 

4 filtering the signal produced by the detector with the variable width digital 

5 filter based upon the determined peak bandwidth. 

1 4. A system for detecting mass to charge ratio of ions, the system 

2 comprising: 

3 a. a mass spectrometer comprising: 

4 i. a vacuum; 

5 ii. an ion-optic assembly adjacent the vacuum; 

6 iii. a free-flight tube adjacent the ion-optic assembly; 

7 and 

8 iv. a detector adjacent the free-flight tube; 

9 b. a control system for receiving a signal from the detector; 

10 and 

1 1 c. means for displaying the signal from the detector; 
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12 wherein the control system includes a variable width digital filter for 

1 3 filtering the signal produced by the detector based upon a desired bandwidth 

14 corresponding to apeak width of the signal corresponding to an ion of interest. 

1 5. The system of claim 4 wherein the control system includes a look- 

2 up table for selecting the desired bandwidth based upon a mass to charge ratio of an ion 

3 of interest. 



1 

2 



6. 

an oscilloscope. 



The system of claim 4 wherein the means for displaying comprises 
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Figure 3 Example spectra data with different digital fillers 
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Figure <-| Detailed view of Peak 1 
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Figure k Detailed view of peak 3 
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